Simulated single molecule microscopy with SMeagol

SMeagol is a software tool to simulate highly realistic microscopy data based on spatial systems biology models, in order to facilitate development, validation, and optimization of advanced analysis methods for live cell single molecule microscopy data. Availability and Implementation: SMeagol runs on Matlab R2014 and later, and uses compiled binaries in C for reaction-diffusion simulations. Documentation, source code, and binaries for recent versions of Mac OS, Windows, and Ubuntu Linux can be downloaded from http://smeagol.sourceforge.net.Comment: v2: 14 pages including supplementary text. Pre-copyedited, author-produced version of an application note published in Bioinformatics following peer review. The version of record, and additional supplementary material is available online at: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btw10

Noise-Induced Min Phenotypes in E. coli

The spatiotemporal oscillations of the Escherichia coli proteins MinD and MinE direct cell division to the region between the chromosomes. Several quantitative models of the Min system have been suggested before, but no one of them accounts for the behavior of all documented mutant phenotypes. We analyzed the stochastic reaction-diffusion kinetics of the Min proteins for several E. coli mutants and compared the results to the corresponding deterministic mean-field description. We found that wild-type (wt) and filamentous (ftsZ( −)) cells are well characterized by the mean-field model, but that a stochastic model is necessary to account for several of the characteristics of the spherical (rodA(−)) and phospathedylethanolamide-deficient (PE(−)) phenotypes. For spherical cells, the mean-field model is bistable, and the system can get trapped in a non-oscillatory state. However, when the intrinsic noise is considered, only the experimentally observed oscillatory behavior remains. The stochastic model also reproduces the change in oscillation directions observed in the spherical phenotype and the occasional gliding of the MinD region along the inner membrane. For the PE(−) mutant, the stochastic model explains the appearance of randomly localized and dense MinD clusters as a nucleation phenomenon, in which the stochastic kinetics at low copy number causes local discharges of the high MinD(ATP) to MinD(ADP) potential. We find that a simple five-reaction model of the Min system can explain all documented Min phenotypes, if stochastic kinetics and three-dimensional diffusion are accounted for. Our results emphasize that local copy number fluctuation may result in phenotypic differences although the total number of molecules of the relevant species is high

Thermodynamic Modeling of Variations in the Rate of RNA Chain Elongation of E. coli rrn Operons

AbstractPrevious electron-microscopic imaging has shown high RNA polymerase occupation densities in the 16S and 23S encoding regions and low occupation densities in the noncoding leader, spacer, and trailer regions of the rRNA (rrn) operons in E. coli. This indicates slower transcript elongation within the coding regions and faster elongation within the noncoding regions of the operon. Inactivation of four of the seven rrn operons increases the transcript initiation frequency at the promoters of the three intact operons and reduces the time for RNA polymerase to traverse the operon. We have used the DNA sequence-dependent standard free energy variation of the transcription complex to model the experimentally observed changes in the elongation rate along the rrnB operon. We also model the stimulation of the average transcription rate over the whole operon by increasing rate of transcript initiation. Monte Carlo simulations, taking into account initiation of transcription, translocation, and backward and forward tracking of RNA polymerase, partially reproduce the observed transcript elongation rate variations along the rrn operon and fully account for the increased average rate in response to increased frequency of transcript initiation

Time-resolved imaging-based CRISPRi screening

Our ability to connect genotypic variation to biologically important phenotypes has been seriously limited by the gap between live-cell microscopy and library-scale genomic engineering. Here, we show how in situ genotyping of a library of strains after time-lapse imaging in a microfluidic device overcomes this problem. We determine how 235 different CRISPR interference knockdowns impact the coordination of the replication and division cycles of Escherichia coli by monitoring the location of replication forks throughout on average >500 cell cycles per knockdown. Subsequent in situ genotyping allows us to map each phenotype distribution to a specific genetic perturbation to determine which genes are important for cell cycle control. The single-cell time-resolved assay allows us to determine the distribution of single-cell growth rates, cell division sizes and replication initiation volumes. The technology presented in this study enables genome-scale screens of most live-cell microscopy assays

Protein Pattern Formation

Protein pattern formation is essential for the spatial organization of many intracellular processes like cell division, flagellum positioning, and chemotaxis. A prominent example of intracellular patterns are the oscillatory pole-to-pole oscillations of Min proteins in \textit{E. coli} whose biological function is to ensure precise cell division. Cell polarization, a prerequisite for processes such as stem cell differentiation and cell polarity in yeast, is also mediated by a diffusion-reaction process. More generally, these functional modules of cells serve as model systems for self-organization, one of the core principles of life. Under which conditions spatio-temporal patterns emerge, and how these patterns are regulated by biochemical and geometrical factors are major aspects of current research. Here we review recent theoretical and experimental advances in the field of intracellular pattern formation, focusing on general design principles and fundamental physical mechanisms.Comment: 17 pages, 14 figures, review articl

ePlant and the 3D Data Display Initiative: Integrative Systems Biology on the World Wide Web

Visualization tools for biological data are often limited in their ability to interactively integrate data at multiple scales. These computational tools are also typically limited by two-dimensional displays and programmatic implementations that require separate configurations for each of the user's computing devices and recompilation for functional expansion. Towards overcoming these limitations we have developed “ePlant” (http://bar.utoronto.ca/eplant) – a suite of open-source world wide web-based tools for the visualization of large-scale data sets from the model organism Arabidopsis thaliana. These tools display data spanning multiple biological scales on interactive three-dimensional models. Currently, ePlant consists of the following modules: a sequence conservation explorer that includes homology relationships and single nucleotide polymorphism data, a protein structure model explorer, a molecular interaction network explorer, a gene product subcellular localization explorer, and a gene expression pattern explorer. The ePlant's protein structure explorer module represents experimentally determined and theoretical structures covering >70% of the Arabidopsis proteome. The ePlant framework is accessed entirely through a web browser, and is therefore platform-independent. It can be applied to any model organism. To facilitate the development of three-dimensional displays of biological data on the world wide web we have established the “3D Data Display Initiative” (http://3ddi.org)

Overview of the MOSAiC expedition - Atmosphere

With the Arctic rapidly changing, the needs to observe, understand, and model the changes are essential. To support these needs, an annual cycle of observations of atmospheric properties, processes, and interactions were made while drifting with the sea ice across the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition from October 2019 to September 2020. An international team designed and implemented the comprehensive program to document and characterize all aspects of the Arctic atmospheric system in unprecedented detail, using a variety of approaches, and across multiple scales. These measurements were coordinated with other observational teams to explore cross-cutting and coupled interactions with the Arctic Ocean, sea ice, and ecosystem through a variety of physical and biogeochemical processes. This overview outlines the breadth and complexity of the atmospheric research program, which was organized into 4 subgroups: atmospheric state, clouds and precipitation, gases and aerosols, and energy budgets. Atmospheric variability over the annual cycle revealed important influences from a persistent large-scale winter circulation pattern, leading to some storms with pressure and winds that were outside the interquartile range of past conditions suggested by long-term reanalysis. Similarly, the MOSAiC location was warmer and wetter in summer than the reanalysis climatology, in part due to its close proximity to the sea ice edge. The comprehensiveness of the observational program for characterizing and analyzing atmospheric phenomena is demonstrated via a winter case study examining air mass transitions and a summer case study examining vertical atmospheric evolution. Overall, the MOSAiC atmospheric program successfully met its objectives and was the most comprehensive atmospheric measurement program to date conducted over the Arctic sea ice. The obtained data will support a broad range of coupled-system scientific research and provide an important foundation for advancing multiscale modeling capabilities in the Arctic

Systembiologisk modellering på molekylär nivå

Implementation and analysis of mathematical models can serve as a powerful tool in understanding how intracellular processes in bacteria affect the bacterial phenotype. In this thesis I have implemented and analysed models of a number of different parts of the bacterium E. coli in order to understand these types of connections. I have also developed new tools for analysis of stochastic reaction-diffusion models. Resistance mutations in the E. coli ribosomes make the bacteria less susceptible to treatment with the antibiotic drug erythromycin compared to bacteria carrying wildtype ribosomes. The effect is dependent on efficient drug efflux pumps. In the absence of pumps for erythromycin, there is no difference in growth between wildtype and drug target resistant bacteria. I present a model explaining this unexpected phenotype, and also give the conditions for its occurrence. Stochastic fluctuations in gene expression in bacteria, such as E. coli, result in stochastic fluctuations in biosynthesis pathways. I have characterised the effect of stochastic fluctuations in the parallel biosynthesis pathways of amino acids. I show how the average protein synthesis rate decreases with an increasing number of fluctuating amino acid production pathways. I further show how the cell can remedy this problem by using sensitive feedback control of transcription, and by optimising its expression levels of amino acid biosynthetic enzymes. The pole-to-pole oscillations of the Min-proteins in E. coli are required for accurate mid-cell division. The phenotype of the Min-oscillations is altered in three different mutants: filamentous cells, round cells and cells with changed membrane lipid composition. I have shown that the wildtype and mutant phenotypes can be explained using a stochastic reaction-diffusion model. In E. coli, the transcription elongation rate on the ribosmal RNA operon increases with increasing transcription initiation rate. In addition, the polymerase density varies along the ribosomal RNA operons. I present a DNA sequence dependent model that explains the transcription elongation rate speed-up, and also the density variation along the ribosomal operons. Both phenomena are explained by the RNA polymerase backtracking on the DNA.Felaktigt tryckt som Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 71